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We consider the vacuum energy of the electromagnetic field in the background of spherically 
symmetric dielectrics, subject to a cut-off frequency in the dispersion relations. The effect of this 
frequency dependent boundary condition between media is described in terms of the incomplete 
(^-functions of the problem. The use of the Debye asymptotic expansion for Bessel functions allows 
to determine the dominant (volume, area, . . . ) terms in the Casimir energy. The application of 
these expressions to the case of a gas bubble immersed in water is discussed, and results consistent 
with Schwinger's proposal about the role the Casimir energy plays in sonoluminescence are found. 
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I. INTRODUCTION 



The Casimir effect ||1|-|3| arises as a distortion of the vacuum energy of quantized fields due to the presence of 
boundaries (or nontrivial topologies) in the quantization domain. This effect, which has a quantum nature associated 
with the zero-point oscillations in the vacuum state, is significant in diverse areas of physics, from statistical physics 
\^ ■ to elementary particle physics and cosmology. 

JC ', In particular, in the last years, there has been a great interest in the Casimir energy of the electromagnetic field in 

LJ the presence of dielectric media, due to Schwinger's suggestion Q that it could play a role in the explanation of the 
f~^ phenomenon of sonoluminescence Iq] . 

[^ — . There are essentially two approaches to the subject: One of them consists in summing up retarded van der Waals 

f^ forces between individual molecules IqM] . The second one makes use of quantum field theory to evaluate the vacuum 
'""' , energy of the electromagnetic field in the background of dielectric media (see e.g. |q-[l^,^,[l3|-[l6[ ) . The relation 
^p ' between these two approaches is not well established. Only for a dilute b all, and up to second order in a perturbative 
(~| \ expansion, both methods have been shown to yield the same answer [ [l7|Jl8| . 

'^ • Regarding the relevance of the Casimir effect to sonoluminescence, the results obtained by different groups through 

O i] several calculation techniques (as Green's functions methods, van der Waals forces, ^-function methods and asymptotic 
O . developments for the density of states - see references p9|- p^ , p7| , p5[ -p^ among others) are rather controversial, and 
"^ ' some basic issues remain to be clarified. 
^ , In particular, there is no agreement about the renormalization necessary to remove the singularities appearing in 

the vacuum energy, a fact that renders the physical interpretation of finite parts difficult. But this inconveniences may 
. , have their origin in the fact that the models usually employed in describing dielectric media mostly do not incorporate 
j^ ' a realistic frequency dependent dispersion relation, then leading to an inadequate ultraviolet behavior. 

In a recent paper ||3l[| , a nonmagnetic dielectric ball with a frequency dependent permittivity (a high frequency 
approximation to the Drude model) has been considered. There, it has been shown that a very simple pole structure 
results for the corresponding ^-function, and only a volume energy counterterm (to be absorbed in the mass density 
of the material) is needed to render the Casimir energy finite. Neither surface nor curvature counterterms are needed. 
With the ultraviolet behavior under control, it makes sense to analyze finite parts of the Casimir energy for realistic 
media. In this context, the analysis of the simple model to be considered in this article is a step in the direction of 
incorporating finite frequency contributions. 

It is our aim to contribute to the understanding of the problem by studying a model which incorporates a frequency 
cut-off fl in the boundary conditions at the separation between dielectrics, to emulate the behavior of real dielectrics. 
To this end, we will assume that this boundary is completely transparent for frequencies greater than Q, for which the 
dielectric constants take the values corresponding to the vacuum. A similar dispersion relation has been considered 
in [3^] for the case of a dilute medium. 

In what follows, we will evaluate the Casimir energy for an arrangement of two such media with spherical symmetry, 
by adding the eigenfrequencies of the system. The existence of a cut-off will allow us to subtract the contribution 
of high frequencies, to > fl. On the other hand, low frequency contributions will be represented in terms of the 
incomplete ^-functions of the problem, as introduced in |p3| . Finally, the asymptotic uniform expansions for Bessel 
functions will allow for the necessary analytic extensions and (for a large cut-off 51) for the identification of volume 
and surface contributions to the Casimir energy. 



These results will be applied to the case of a gas bubble in water, a situation of interest to sonoluminescence. 
As we will see, our results seem to support Schwinger's proposal about the role played by Casiniir energy in this 
phenomenon. 

The paper is organized as follows. The model is presented in Section [[|, where the analytic extensions of the incom- 
plete ^-functions are constructed. The (finite) number of modes effectively contributing for each angular momentum 
is evaluated in Section III. The bulk and surface contribution to the Casimir energy are obtained in Section IV. In 
Section M we evaluate the electromagnetic pressure on the separation between dielectrics and the change of vacuum 
energy with respect to the volume. The application of these results to sonoluminescence is also discussed there. In 
Section W^ we present a summary and discussion of our results. Finally, in Appendix H the representation of the 
vacuum energy as an integral on the complex plane which we are employing is justified from a mathematical point of 
view. 



II. THE MODEL AND ITS INCOMPLETE C-FUNCTION 
A. The model 



Our aim is to evaluate the Casimir energy of a spherical dielectric ball or bubble of radius a and indices (relative to 
the vacuum) /ii(w), ei(w), immersed in a second medium of indices /^2(i^), e2('j-')- As a rough model of real dielectric 
media, we will suppose that /ii(w) and ei(w), with i = 1,2, are constants up to a common cut-off frequency Q, 
while their values for lu > H, are those of the vacuum. This last restriction reflects itself in a frequency-dependent 
boundary condition for the electromagnetic field at the separation between the dielectric media, making the boundary 
completely transparent for those modes of frequency greater than the cut-off f2. 

We will evaluate the Casimir energy of the electromagnetic field in this arrangement of media by summing its 
eigenfrequencies up to the cut-off, i.e. over all uj < il. Notice that disregarding, under the above mentioned conditions, 
the contribution of modes with frequency lu > Q amounts to redefining of the zero energy level by the subtraction of 
a (divergent but) a-indcpcndcnt quantity. 

Consequently, for oj < fl the electromagnetic field satisfies dielectric boundary conditions at the surface of the ball, 
which in spherical coordinates lead to 



Es,'p\r=a+ - -^e,0lr=a- 



1 

Ail 
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/^2 



Inside the dielectrics, the electric field satisfies the Helmholtz equation. 



(1) 



(2) 



and similarly for the magnetic field, B . One can consider the transversal electric (TE) modes, taking the electric field 



as 



^;,m = //(r):^y;,™(0,0), 
and separately the transversal magnetic modes (TM) , with the magnetic field given by 

with ^ = 1, 2, ... in both cases. In the previous equations 

L = —i~r^ X V = —ifpde + i9— — -d^. 
For the TE modes, the boundary conditions (1^) imply 



sm( 



fi{r)l=a+ ^ fi{r)\r=a- ^ - ^r [r/i (r)] 



- —dr[rfi{r)] 



(3) 



(4) 



(5) 



(6) 



For the TM modes, the same conditions reduce to 



— giir) 



= —9i{r) 



dr[rgi{r)] 



.+ ^^l £1 



dr[rgi{r)] 



(7) 



Then, for example, we get for fi{r) 

--^ r/u^)] ^:5 — fi{n = -Mi.2ei,2-^//(r), (8) 

for r ^ a which, together with the boundary conditions in Eq. (|g), imphes that fi{r) is a continuous piecewise 
differentiable function with a discontinuous first derivative at r = a. 

In order to have a discrete spectrum, we enclose the system inside a large concentric conducting sphere of radius 
R ^ a, obtaining also the Dirichlet condition, fi{r) = at r = R, for the functions in the domain of the relevant 
differential operator. We will take the R -^ oo limit at the end of the calculation. 

In Appendix H, we show that the eigenfrequencies corresponding to TE modes are determined by the zeroes of the 
function 



ATE 

^1+ 



l/2{z) = Jl+l/2{zi) [yi+l/2{zo)Ji'+T^/2{z2) " Ji+1/2 (zo)X+i/2 (^2) | 
C'^/+l/2(^l) {yi + l/2izQ)Jl+l/2iz2) - Jl + l/2iza)yi + l/2{z2)} , 



(9) 



where 



Ji+i/2{w) = w ji{w) = J—Ji+i/2{w)., yi+i/2{w) = w yi{w) = J —yi+i/2{w) (10) 



are the Riccati - Bessel functions. In Eq. (|9|), z = a{uj/c), z\^2 = z^t\;i[i\;i^ zq = zR/a^e2fJ-2, and ^ = ^/^i^. In the 

same Appendix we show that the zeroes of the function Af_^^,2{z) in the open right half plane of the variable z are 
all real and simple. 

For the TM modes, the same analysis can be done for the function gi{r) = — gi{r), defined for r ^ a and satisfying 
the Neumann boundary condition, g'lir) = at r = i?. In this case, the eigcnfrcquecies are given by the zeroes of the 
function 

Ar+f/2(^) = ^^ + 1/2(^1) {X+l/2(^2)^/+l/2(^0) - J'/+l/2(^2)X+l/2(^0)} 

1 ^^^^ 

-jJUiI2{z\) |3^;+i/2(^2)J'/+i/2(^o) - J;+i/2(^2)X+i/2(^o)} 

contained in the open right half plane of the variable 2, which are also real and simple. 

In order to simplify our calculations, in what follows we will consider both media to be nonmagnetic (/ii = 1 = /i2 
for all frequencies), while keeping y^ = n\ and ^JT^ — 712 arbitrary for a; < fi. 

B. Vacuum energy and incomplete (^-functions 

We will be interested in evaluating differences between vacuum energies corresponding to situations differing in the 
value of the radius a. Then, as remarked above, we can disregard the contributions of those modes with frequencies 
w > J7 since, being independent of the position of the boundary (and also of the low frequency refraction indices), they 
cancel out (whatever the regularization employed in defining the vacuum energy would be). This simply amounts to 
performing an a-independent subtraction, which is nothing but redefining the zero energy leveln. 

Therefore, for the TE modes, we must evaluate the (finite) sum 



^For example, in the framework of the (^-function regularization, the contribution of the frequencies greater than the cut-off 
can be defined as the analytic extension to s = —1 of a series convergent for JR(s) large enough: 



E^^{a)=^ ^ 21/ ^ -hio^^n^— J2 ^ ^ ^'^,n, (13) 

i/=3/2 n=l i/=3/2 n=l 

and a similar expression for the TM modes. In Eq. (13), N^ is the number of positive zeroes of A^^(z), z^.„, which 
are less than or equal to a; = aVi/c. The degeneracy due to the spherical symmetry is 2v — 2/ + 1, and vq is the 
maximum value of v for which Ny > 1 . 



We are interested in an analytic (rather than numeric) evaluation of Eq. (13). So, although this is a finite sum, we 
will employ the summation method developed in [ p3| , based on the evaluation of an incomplete C- function. We can 
use the following representation: 



n— 1 n—1 



(14) 



s=-l 



where the sum in the right hand side exists as an analytic function^of s G C. 

Since AT^(z) has only real zeros in the open right half z-plane, and its non- vanishing zeros are all simple (see 
Appendix |^, we can employ the Cauchy theorem to represent the sum in the r.h.s. of ( p^ ) as an integral on the 
complex plane. 



N^ 



n=l 



.-.A^^'(z) 



E a = ijz: f ^-'T^fEhi ^^' (15) 



2m Jc AV{z) 



where the curve C encircles the first N^ positive zeros of A^^(z) counterclockwise. 

For 3fJ(s) large enough, the contour C can be deformed into two straight vertical lines, one crossing the horizontal 
axis at ^{z) = x and the other at 5R(z) = 0+. Indeed, the integrand can be expressed in terms of modified Bessel 
functions through the substitutions p4| 



Jy (e'f w) = e'f "^ h (w) , Y^ie'i w) = e*2 (^^+1) /,, [w) - | e-'i" K^ {w) , 
relations valid for — tt < arg(fi;) < 7r/2. Thus, we get (see Eq. (pi)) 



(16) 



,TE/ . N o (^^+1/2) jniTTR 3 {^^,n2Rv, 



Al^{x + iy)^-e2 ^ ' ^ ^^ „i ^ ;^^(^_11) [_„^ j^(„2 «) /^(„^ „) 
+n2 Iu{niv) Il(n2v)] +/,,( ) [rii Ky{n2v) ll{niv) - n2 Iy{niv) Kl{n2v)] 



(17) 



where v — {—ix+y). Taking into account the asymptotic behavior of the modified Bessel functions for large arguments 
P4| , it is easily seen that, for < x 7^ Ziy,„,Vn, the integral 



1 ^2;+joo \TE'{^\ 



C^{s,x)^— z-^-^p^dz, (18) 



Eo{s) = lhnJ22,. ^ ^) , (12) 

where u = I + 1/2, and ujI^„ are the zeroes of A'^^ (auj/c) and A'^'^' (auj/c) taken with ni — 1 — n2, i.e. in a situation 
indistinguishable from one where the external sphere contains only vacuum. Obviously, Eo{s) is independent of a, ni, and n2. 
It can only depend on R and Q.. 

^Notice that the sum in the r.h.s. of Eq. ([L4|) evaluated at s = gives A^^, the number of eigenfrequencies contributing 
effectively to the Casimir energy of the field for a given value of the angular momentum I = v — 1/2, once the a-independent 
subtraction adopted to define it has been made. 



converges absolutely and uniformly to an analytic function in the open half-plane 5R(s) > 1. Without loss of generalit 
and for calculational convenience, we will restrict ourselves to real values of s, and evaluate the function in Eq. (|l 
on the half-line s > 1, from which it can be meromorphically extended to the whole complex s-plane. 
Therefore, for s > 1, 

E ^.-: = cr(^,o+)-cr(s,x). (19) 

n=l 

Moreover, since the left hand side of ( |l9| ) is holomorphic in s, the singularities of Cj-^(s,a;) must be independent of 
X. In particular, this allows us to write the vacuum energy as the analytic extension 



Tic 



n ^ — ^ ■- 



(20) 



,.. , , . ... (s,x} 
a ■^ 

Entirely similar conclusions arc obtained for the TM case. 
Moreover, for ?R.{w) > we have Q| 

(21) 
II {e-'-w) = e-^^^"" - 1) {Il{w*))* , Kl (e"'" w) = e^^i"" + 1) {Kl{w*))* + ivre*'^ (/^(w*))* • 



So, changing the integration variable in Eq. ( |18[ ) by z -^ iv^^) x, and calling t = {y~ i)x/v, we can straightforwardly 
write for real s > 1 



.T^.„.^_«J- -*?(^ + l) r"'%-.rf(lnA^^(*^t)) 



Cl^(s,x)=^l^e 2^ 'J t-^ ^_^^rfM, (22) 

where we have now called z = xjv > 0. 

A similar expression is obtained for Cj*^(s, x), corresponding to the TM modes. 

C. The analytic extension of incomplete (^-functions 

In order to construct the analytic extension of Cj^('S, x) to s ~ — 1, we will subtract and add to the integrand in (p2) 
the first few terms obtained from the uniform asymptotic (Debye) expansion [Q of the Bessel functions appearing in 
the expression of Aj^^{ii't) (see Eq. (O)), which is valid for large ly with fixed t: 



-^ ^-=D^'^{t) + 0{,y-^), (23) 

where 

D^^it) = ^i?W (t) + D^^lit) + iy-'D^r^'\t). (24) 

In the above expression, the functions Drp'^[t), k = 1,0,-1, explicitly shown in Appendix H, are algebraic functions 
of i. Notice that we have discarded contributions coming from terms containing K^{ ^'^ Ri^t -^ ^^ ^^^ derivative (see Eq. 
(|l7|)), since they vanish exponentially when i? — > cx). We will see that this approximation allows for the identification 
of the volume, surface, . . . contributions to the vacuum energy. 
So, we must consider the integral 



.oo-« d(lnAl'^(iiyt)) f°°-'^ „„ f°°-'^ fd(lnAP(iiyi)) 
/ t-'^ V^ '-^dt= t-'Dl^{t)dt+ t-' {^ v^ ^ 



Dt^it)^dt. (25) 

The second integral in the right hand side of Eq. (|2^) converges for s > — 2, since the integrand can be estimated 
by means of the contribution of the next (©(i/^^)) term in the Debye expansion (see Eq. (E3)), which behaves as 
©(t^^^"^) for large \t\. This term could be evaluated numerically at s = —1. This will not be done in this paper. 



The Debye expansion applied to the TM modes case gives 

Dl''\t) = z.4^],(i) + 4"],(t) + :^-'D^^-^^it), (26) 

instead of Eq. (p4) , and leads to a decomposition similar to Eq. (E3) . The algebraic functions D^lj (i) are also shown 
in the AppendixlH. 

In what follows, we will evaluate only the first integral in the right hand side of ( P5[) (and the analogous expression 
for the TM modes), retaining only those terms of its expansion in powers of i/^^ which are consistent with the 
approximation made in Eq. ( P3h . 

Notice that the integrand, IJ^^{t) (D^^^ (t) for TM modes), is an algebraic function behaving as 0{t'^) for large |i| 
(see Appendix ^). So, the integral converges absolutely and uniformly for s > 1, where it defines an analytic function 
which can be meromophically extended to the region of interest of the parameter s. As we will see, this extension 
reveals the singularities of (^'^^{s,x) (C^^{s,x)) as simple poles with x- independent residues (a necessary condition 
to get a finite result for any s in Eq. (E9[)). Notice that this statement must be valid for the contribution of each order 
in V. 

We begin the calculation by considering th e te rms of dominant order in the Debye expansion^, v Dj,^{t). By virtue 
of the analyticity of the integrand (see cqs. (|Bl|)), for s > 1 we can deform the path of integration to write 



r'vD^^l{t)dt^v t-^-1 Vl + m^t2-v/l + n2^^^ + Ji+ ' , dt 



na^ i?2 ^2 



-/ 



1 1 



Cn J^ 1 L g \ 
"2 -ft , m 712 ^ 112 fl 1 J. 
n-i - "2 H V — 7r3 ^— dt. 
a IV j 

The first integral in the r.h.s. of eq. (p7|), which contains the whole dependence on x = j/z, is holomorphic in s and 
can be directly evaluated at the required value of this parameter. On the half-line (1,cxd) we have subtracted and 
added the first terms in the series expansion of Drp^(t) for large t, which makes the second integral convergent for 
s > —2. The third one must be evaluated for s > 1 and then analytically continued to the relevant values of s. This 
can be exactly done, and its contribution to C,^^{s, x) in Eq. (E3) is 

j,^ I i/^^* cos (I (1 + s)) I 2nin2R {aril - an2 + n2R) ~ (a^ ni) + ani R ~ an2 R\ 

^iCi/ (s, 2;) o: = - — - I - ; • (28) 

This expression is analytic at s = and has simple poles at s = ±1, which are the only singularities of the dominant 
contribution to Cj (s, a;) in this asymptotic expansion for ^{s) > —2. In particular, its residue at s = — 1 is 

n A /-TEr M ly'^ {ani + {n2-ni) R) 

ResAiC^^(s,a;)|s=-i = . (29) 

2 ni 77,2 TTil 

Notice that, as anticipated, up to this order in the Debye expansion the residue is independent of x. 

For example, for the dominant contribution to (J^{s,x) (which coincides up to this order with Cj*^(s,a;) - see 
footnote^), 

cr(s,x) = Aicr(s,x) (i+o(t.-i)), (30) 

one straightforwardly obtains the Laurent expansion around s — ~1 



^ Since Dj,j^^{t) is identical to D^l,{t) (see eqs. (Bl) and (B5)), we obtain the same result for the dominant contribution to 

C™(s,x). 



2 ni ^2 TTit (1 + sj 

f 2 

-5R < {aril + R (n2 — ni) (1 — 2 ni ^2) +2ni n2 R (119 R + ta z ( yl + e^'^^rii^ z^ 

I4nin2 7ri? V" V 



e '■^n2^i?^2:^ 



-v/l + e-*'^n22 2:2 + A/lH ^5 +2 {ani log{n2 R) - ni R log{n2) 

+n2R log(ni) )log(2)+ (71,2 - ni) (2 711/12) R- 2nin2^ 2 (arii + (-ni +712) i?) log(j/) 

a 

—2 77-2 -R log (~iniz + v 1 + e^' ^^rii^ z^ J +2 77,1 _R log I — i r72 z + yl + e^' ■^n22 z^ 1 



(31) 



-2 a 711 log I ~in2Rz + a^l + ^ I +0(1 + 5). 

On the other hand, for a; — > 0+ a similar calculation leads to 

A,cr(.,. . 0+) = Aicr (5,x = 0+) ^ "' ^""^ + ("^ - "^) ^) 



4ni 7T,2 TI" 



2r7,l 772 TT _R (1 + s) 

(2 711 log(2n2)-2?i2 log(2ni) + (7ii-7i2) (1-2 log(7/))) + 0(5 + 1). 



(32) 



Similar calculations are required in order to get the contributions to C,, (x, s) ( C,, (x, s)) coming from the next 



to leading terms in Eq. ( p4| ) (Eq. ( |26| ) respectively). In fact, from Eqs. ( |B2[ ) (Eq. ( |B6[ ) for the TM case) it is easily 



seen that Dj,'^{t) {D^lj{t)) ^ t ^. So, its contributions to the first integral in the right hand side of Eq. (M) (or 
the equivalent for TM) converges to an analytic function for s > —2, and does not affect the residue of Cj^a;, s) 
{C™{x,s))a.ts = -1. 

For the second order term in the Debye expansion of the difference [Cj^(sj x = 0"*") — C,J^{s, x)\ around s — —1, a 
straightforward calculation leads to 

AoCj^(s, 0+) - AoCj^(s, x) = = u (--^ e(ni x - 7.) - -^ 8(712 x - z/) - -^ e{n2Rx/a - v) 

( 4ni 4 712 4 712^1 

(33) 



2 „ . ^ „ I (n-^ + ri< ) (;/> — u) ii> — 7i^ \ 
+ — ^e(7.>-7.)i^ arcsinWp i)\_ . \_ +0(s + l). 

7r(71>+71<) \ y (71> - 71<)(7/> +U) ?1>+71</J 

Here, Q{w) is a step function (vanishing for w < and equal to 1 for w > 0), F{(p, k) is the elliptic integral of the 
first kind, 7i< (7i>) is the min(max) {711,712}, and u — maxJT^, i^<}, with u^ = n^x and 7/> = 7i>a::. 
Similarly, for the TM modes we get 



AoCj*'(s,0+) - AoC^Hs,x) = 7. <j --^ 9(711 a; - 7.) - -^ 6(112 x - u) + -^ 6(712 i?x/a - '^) 

4711 4772 4 712 it 



I/) n 8 



+ _:<_e(7i>x-7/)nlarcsiny-|— ^,l--|,^l--|j|+0(5 + l). 



(34) 



where n((/3, 77, /c) is the elliptic integral of the third kind and u = maxjz^, 7^<}. 

Finally, notice that Dj,^ {t) and Drp^J{t) w i^^ (see Eqs. (B3) and (p37|)). So, they do contribute to the poles at 
s = —1. We get 

A_iCj^(5, x) = A_ie*^(5, x)= , "^",7^ , - ^-— - + 0(5 + 1)°, (35) 

8nin27r(l + s) 8772 7ri? (1 + s) 



although the finite parts (which will be not given explicitly here) are different for TE and TM. Notice that the residues 
are independent of x, and cancel out when the difference in Eq. (^9|) (or the equivalent for the TM case) is taken. 

In the followin g Se ction we will evaluate N^, (i.e., the number of modes contributing in Eq. (|l^)) as a function of 
v, and in Section [V we will get their contributions to the vacuum energy. 



III. THE NUMBER OF CONTRIBUTING MODES 



In this Section we address ourselves to the determination of Pq in Eq. (|13|), i.e., the maximum value of v for which 
Ni, > 1. As before, we will follow the method established in p3^ . 
First, notice that 



N^ 



N.{x) ^ Y. 



n=l 



- l/-TEf^ n+ 



[cr(«,o+)-cr(.,x)]|^^„ 



(36) 



s=0 



is a step function of x, having a discontinuity of height 1 at each positive zero Zi/^„ of the function A^^(z) in Eq. 
Then, 1^0(2;) can be determined from the condition 



N^.ix) ^ N^,iz^,^i + 0) = 1, 



(37) 



with N^„{z^„^i - 0) =0. 

For the dominant order of D^^{t), Eq. (Bl), taking into account Eq. ( |2q ) and the fact that the second integral in 
the r.h.s. of Eq. ( |27[) is (finite and) real at s = 0, it is straightforward to obtain, from Eqs. (|2^ ) and (p7|), that 



AiCj^(s = 0,.x) = -^-5R{^(yr 



e *^ni^z^ ^yl 



e *^n2^z^ 



1 + e-"^ 



n2 U z^ 



/ ni'^ B? z"^ 

log(l + Vl + e-*^ ni2 z2) + log(l + Vl + e-*'^ na^ z2) -log(l + Wl + e^*'^ ^ ) 



(38) 



In particular, for x ^ 0^ 



AiCr(s = 0,x = 0+) 



Similarly, from Eq. (B2), a straightforward calculation leads to 

AoCP(s - 0,x = 0+) - AoCP(,s = 0,x) = 

= 5R I— nog(l + e-'^'ni^ z^) + log(l + er'^'n^^ z^) + log(a2 + er'^'n-,^ R^ z^)] 
|_ 47r 

Now, calling 

iVj^(x) ^ Aie^(s = 0, 0+) - AiCj^(s - 0, a;) + AoCj^(s = 0, 0+) - AoCj^(s = 0, a;) 

= iV,(x)+ 0(^.-1), 



(39) 



(40) 



(41) 



we see that it gives a smooth approximation to the step function in ( p^ ) for v ~:^ 1. So, following p^, we will 
approximate 1^0(2;) by the value of v for which N'^^{x) = 1/2. 

Since z = x/v (with x = af2/c), it can be easily seen that Nj^{x) = for i' > n2Rx/a (=^ 1 > n2Rz/a > ^1^2-^), 
while for n2Rx/a > v > ni_2X we have 



Nr{x)=- 



n2^R^z^ / /n22i?2^2 
;; 1 — arctan a / 1 



(42) 



Now, we will assume that 7Vj-^(x) = 1/2 for vq < n2Rx/a, write e 



2 _ n2 R a: 



1, and determine e iteratively 



from the series expansion of the right hand side of Eq. (|42|) around e = 0. This leads to 



rp^ 712 i? n J ,TE ( ^'2R^\ , iTE ( '"''^R^\ I /o /^ -^ ^ 



•'' -^1^ ['-''- y-^} ^'i-y-^} +''1—1 ■>• («' 

where 

, TiS _ ol/3 :i^^ i,TB _ o2/3 '^ ^ ^44^ 

^=1 -^ 4 21/3' ^2 -^ 160 22/3- ^^^^ 

Notice that this result does not depend on a or ni, i.e., the radius of the bubble and its refraction index. 
For the TM modes, a similar calculation shows that N]^^ [x) = N'^^{x) + 1/2, and 

^TM^]^L_^TM hl^Y'\^TM ^-^^^^^ +0\ — ] ), (45) 

with 

" -^ 4 21/3' ^2 -^ 160 22/3- ^^^> 




,TM 



is also independent of a and ni 



IV. THE DOMINANT CONTRIBUTIONS TO THE VACUUM ENERGY 

In this Section we will evaluate the contribution to the vacuum energy due to the dominant orders in the Debye 
expansion of incomplete C- functions, obtained in Section n. 

A. Bulk contributions 

According to the results in Section ||, we need the Laurent expansion of Ai(^J-^(s, x) around s = — 1 (given in Eqs. 
( pi| ) for arbitrary X, and in Eq. (B3) for a; = 0+). As already remarked, the contribution of singular parts to the right 
hand side of Eq. (n9| ) cancel out, since the residues are independent of x (see Eq. (p9[)). For the difference of the finite 
parts for v < Vq'^ we get 



[AiCr(.,0+)-AiCr(.,x)]|^_^ = 
1/2 z {Q{nix - u) G{niz) - e{n2X - v) G{n2z) +Q{v'^^ - v) G{n2Rz/a)] 



(47) 



G(^.) = ^^^IH_M!fi±^ZZi). (48) 

27r 271-01 



where 



Replacing z — xjv and x — aQ,/c^ we obtain for the dominant order contribution to the sum in Eq. (P0[), 

i^<niaO/c i/<n2ar2/c u<uT^ 

where v = I -\- 1/2, with / = 1, 2, . . . Since ^'^G{1/ ^) has a pronounced maximum at /^ « 1/2, the approximation to 
E^^{a) given by Ai£^-^^(a) is justified as long as aVl/c » 1. 

Notice that the last term in the right hand side of Eq. ([49|), the only piece which is a function of i?, depends neither 
on a (the radius of the ball) nor on ni (the refraction index of the internal medium). In fact, as shown in the previous 
section, v^^ is independent of these parameters (see Eq. ([4^)). Therefore, this is the only piece remaining in the limit 
a — > (no internal bubble), where the other two terms are vanishing. 

In order to sum up the first and second terms in the right hand side of Eq. (B^) notice that, even though fk{v) = 
v'^G{nkaVt/vc), A; = 1 or 2, and its first derivative are finite and bounded for v g [IT^^fc] (with Vk = [n^aO/c— l/2] + l/2, 
where the square bracket denotes the integer part), the second derivative is unbounded near Vk- This is so because 



/fcM-gfeH + OK-i^)^ (50) 

where 

g,(^) = i^K-^)i. (51) 

Therefore, we can subtract and add gk{v) to fk{i^), and apply the Euler - Maclaurin summatfon formula [ p5[ to the 
difference, 



J2 [/feM ~ 5fc(^)] = / ifkix) - 9kix)] dx+- ([/fc(3/2) - g,.(3/2)] + [AK) - g^iiyk)]) + 

.~r,^ J3/2 ^ 



^=5/2 -^"^-l 



Since the second derivative in the argument of the last integral is non-positive, it is easy to see that the remainder 
(this last term) is 0{nkafl/c). Then, a straightforward calculation leads to 

^ rr , , , M (5-16V2) frikanY , 2^/2 f nkaCtY , .^ f rikail 

1^=3/; 

On the other hand, 



i/=3/2 ^ ^ ^ ^ ^ '^ 



1^=3/2 ^ ^ ^ ^ ^ 



(54) 



where a/c = (nfcai7/c — 1/2)— [rikafl/c — 1/2] G [0, 1), and the asymptotic expansion of the Hurwitz ^-function, C(s, v), 
for large v |]35| has been used in the last step. 

Finally, adding the results in Eqs. (m) and (M) (notice that the surface contributions cancel out), taking the 
difference for k = 1,2, and adding a similar expression coming from the third term in the right hand side of Eq. (|4S|), 
we get 



nn ' ' 12 TT \ c / \ c 

(55) 



The same result, with /cff -^ fcff^, is found for Ai_E^*^(a). 

Eq. (pq) shows that the dominant contributions to the vacuum energy in this asymptotic expansion are volume 
terms, in agreement with the claim in P, p3| , p7|j2^ . 

There is a term proportional to the volume of the accessible space (~ i?^), with corrections depending on fractional 
powers of R induced by the cut-off imposed (see Eqs. (^) and (|45|)). These corrections are independent of ni and 
a. There is also a bulk contribution proportional to the volume of the ball (~ a^), which is twice the one obtained 
for the scalar field case discussed in |33|, multiplied by (ni^ — n2^)- So, it is the sign of the difference (nf — 7^3) 
which determines the vacuum energy behavior with respect to the radius of the bubble. In particular, it vanishes for 
rii =712. This will be further discussed in Section M. 
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B. First finite-size corrections 



In order to incorporate the first finite-size correction to tlie vacuum energy, we need tlie Laurent expansions of the 
next to leading order in the asymptotic expansions of the C- functions around s = — 1, quoted in Eqs. (pSh and (p4). 
For the TE case we get 



1 
Ml 



AoE'^ia 



afl 



Y: z^[AoCr(s,0+)-AoCr(s, 



1^=3/2 



c 



4n<- -^^ 4n> ^-^ ^noR ^-^ 



^'<Z-'< 



f <Z.'> 



v<vl 



(56) 



7r(n> + n<) 



K 



n> + n< 



E 



y^ v^F arcsin . 



V^ <1'<U 



'(n> + n<)(i^> — i/) riy — n^ 
(n> — n<)(i/> + v) ' n> + 7i< 



where i/ = / + 1/2, with ^ = 1, 2, . . ., z/< = n<x, i^> = n>a; and K{k) — _F(7r/2, A:) is the complete elliptic integral. 
Similarly, for the TM modes we have 



1 

Ml 



AoE' ^ (a 



aU, 



^ ..[AoCj*^(.,0+)-AoCr(.,x)] 



1^=3/2 






f <!-'< 



1/<I.'> 



v<vl 



nif-:iV'-tjj:^^ 



(57) 



+ 



E M 



I/< <I^<Z^> 






For simplicity, let us assume that the refraction indices are such that (i^< — 1/2) and {y^ — 1/2) are both integers. 
This does not lead to any loss of generality in the result we are looking for since, as in the case of the bulk contributions 
previously worked out, the fractional parts a^ = UkX — [ukx] have no effects on the leading terms of the sums for 
{an/c:$> 1). 

Almost all the sums appearing in the right hand side of Eqs. (pq) and (m) can be trivially solved, since 



^f 



E-^ 

i/=3/2 



+ 0(^/)- 



(58) 



The exception are the last terms appearing in those equations. Once again, these contributions can be approximated 
by means of the Euler - Maclaurin summation formula. 

In so doing, one should remark that the functions in the argument of these sums, say /(j^), vanish (in both cases) 
as a square root at z/< and at Vy. It is sufhcient for our purposes to subtract a function g{h') behaving the same way, 
in order to obtain a difference with a bounded positive first derivative, to which we can apply the Euler - Maclaurin 
formula |3^], 



E (/H-.9(^))= r (/(^ 



) - giu)) dy+ 



(59) 



+ (/(z.< + 1) - g{v^ + 1)) + r {x - [x]) (fix) - g'{x)) dx. 



One can easily verify that the remainder (the last term in the right hand side) is ©(ail/c)^. On the other hand, the 
sum X^jy^jy +1 9{^) can be solved in terms of Hurwitz (^-functions. 
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Gathering these results together we get 



for the TE modes, and 



1 A j,TM.._ 1 (">'-«<') fan\\^^fan 



nn ^ ' 12 n->2 +n^2 I „ ; I „ 



> 



c / \ c 



^(=f^)"l--™(=^)""l-(^ 



(61) 



for the TM case. 

Notice that for the TE modes there are no surface contributions coming from the interphase between dielectrics, in 
agreement with the results in |g2[ for nonmagnetic media. On the other hand, there is a negative surface contribution 
from the TM modes, vanishing for ni = 712- Moreover, in both the TE and TM cases, there are surface contributions 
corresponding to the external boundary (~ R^), which differ in sign and cancel out when added. This is in accordance 
with the fact that Dirichlet and Neumann boundary conditions contribute with opposite surface terms |£^ (see Section 
P). Finally, there are also corrections which depend on fractional powers of R, which are induced by the cut-off imposed 
(see Eqs. ^^ and (|4|)). 

Notice that the procedure followed to evaluate the Casimir energy can be continued up to any given order in the 
asymptotic expansion in Eq. (p^), to get the result up to the corresponding order in powers of (afi/c)"^. 

V. THE CASIMIR ENERGY 

A. The electromagnetic vacuum pressure on the bubble 

Let us stress again that the contributions depending on R in the right hand side of Eqs. (p5h (and the corresponding 
result for the TM modes), (|6^) and ( pi] ) are independent of a and rii. They are exactly canceled if one refers the 
energy to that of the medium with 712 filling completely the interior of the external sphere, by subtracting the same 
expressions with a = 0. In this way, any reference to the exterior radius R disappears and one obtains 

J_ (m^ - 712^) fan\ ^ _ l_ {ni^ -n2^f f an\ ^ fan 

nn^^"'' QiT \c) 12 ni2 + n22 \c)^^\c 



The second term in the right hand side of Eq. (p3) (a surface contribution, qualitatively similar to the one obtained 
in ||32| ) is negative, while the behavior of the first one (a volume term) depends on the sign of [ni^ — n2^). For a 
large cut-off, the volume contribution is dominant, in agreement with the claim in p|, |23| , p7| , p8[ . 

As expected, £{a) ~ for ni — 712. If, for example, 712 > rii, then £{a) is a negative function, monotonically 



decreasing with a. But, as remarked in |29|, the values of £{a) for different ball radius refer to configurations with 



different amounts of material media, and are not directly comparable. 

Instead, we will retain the whole dependence of the vacuum energy with R, as obtained in Eqs. (|5^), ( pOl ) and (|6l|), 
in order to allow for a variation of the refraction indices with the volume of the bubble, while keeping the number of 
molecules of each dielectric constant. This condition is equivalent to demanding that [^ 

[ni(a)2 — 1] a^ = constant, \n2{aY — l] {R^ ~ a^) — constant, (63) 

which implies that 



("i(a)^-l) , 3a2 ("2(0)^-1) 



^^ ' 2a ni(a) ^^ ' 1 {R^ - a^) n2{a) 
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These derivatives, replaced in the expression of the pressure acting on the boundary between dielectrics due to the 
electromagnetic field, 



P{a 
-1 d 



-1 



9^ + "i("W + "^^"W 



Ecas. (a) 



4 TT a^ da 
straightforwardly lead to 



(AiE^^ia) + A,E™{a) + AoE^^{a) + A„E™{a) 



(65) 



p{a) fn 
nn \7 



+ 



1 

^ Tfi~~2 { n2{af-ni{af -3 n2{a) - ni{a) | 
3 (ni(a) — Ij (n2{a) — ni{a) j (ni{a) +8712(0) 



\ni{c 



"-2(a)' 



24 TT 



}+0 



a 17 
c 



c 



fni(a)^ - n2(a)^ 
7ii(a) +712(0)^ 



(66) 



O 



i?Sl 



-2/3 



where the limit i? — > 00 can be safely taken. 

Notice that, even though 712(0) has a tiny derivative (~ R^^), it enters in a term with a large coefficient (~ i?^), thus 
giving a finite contribution to P{a). On the other hand, those terms containing lower powers of R do not contribute in 
the limit R — > cx). Therefore, the expression obtained for P{a) is presumably independent of the boundary conditions 
imposed on the field at r = _R. 

The first term in the right hand side of Eq. (|66|), coming from the bulk contribution to the vacuum energy (Eq. 
(pq)), is clearly dominant for afl/c >> 1, while the second term, coming from the surface contributions (Eq. (61)) is 
less significant in this region. 

The pressure P{a) behaves in the following way, depending on the values of the refraction indices: As expected, it 
vanishes for 712 — ni, and its derivative with respect to the exterior index is negative. 



dPja) 
dno 



16 7r2 



{nl l) 



of^ 



<0, 



(67) 



for 712 > 1. Therefore, for 772 > 7ii, P{a) < and the bubble tends to shrink, while for 712 < 77i, P{a) > and the 
ball tends to expand. In this way, we arrive at the nice picture of a dielectric tending to fill empty space, even for a 
vanishing electric field. 



B. The Casimir Energy 

The Casimir energy as a function of the bubble's radius (for given amounts of dielectric materials) can be obtained 
by integrating —P{a) with respect to the bubble's volume. As remarked above, 772(a) is essentially constant, since 
this dielectric has a very large available volume (see Eq. (|64|)). On the other hand, when the bubble originally filled 
up with a dielectric of index 771 is expanded from a volume Vb = 47raQ/3 to a volume V = 47ra^/3, one finally gets a 
refraction index given by (see Eq. (|63|)) 



7ii(a) = Wf + (772-1) 



v 



Retaining only the dominant term in the expression of the pressure, Eq. (Em , we get 



Ecas. (a) - Ecas. (oq) = - / P{a} dV = 

JVo 



c 



16 7r2 



(?^2 



37^2) {V ~ Vo)+ +g^ 



^(1 + (V-1)^ 



m' Vo 



(68) 




(69) 
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Notice that, up to this first order, the dependence on rii and 712 appears in separate terms. 

This resuh is more easily analyzed in the case of a bubble containing a dilute medium, i.e., when nf — l = ei — l<Cl. 
In this case, we get 

EcasXa) - EcsXao) - nn (^)' | ^^^ ' fj + "^^ (V - Vo) - ^^^^ ^ (V - V,) + 0{nl -if^. (70) 

If the first term is dominant (with ^2 > 1) the Casimir energy increases with the volume. But if, for example, the 
ball's exterior contains just vacuum (^2 = 1), the first term in the right hand side of Eq. (^0|) vanishes, the second 
(with a negative coefficient quadratic in (ei — 1)) is the dominant one, and the Casimir energy decreases with the 
volume as {Vq/V — 1). 

C. Application to Sonoluminescence 

One reason for the great attention devoted to the study of Casimir energies of dielectric media with spherical 
symmetry is the suggestion made by J. Schwinger B that the Casimir effect might have a key rol in the explanation 
of the peculiar phenomenon of sonoluminescence (consisting in the transduction of sound into light, see [^ and 
references therein). 

This phenomenon (for a review, see |g|l) is characterized by the fact that the energy enters a fluid as a sound wave 
(of ~ 26 kHz) which induces on a single gas bubble (air with some proportion of a nobel gas), trapped in a velocity 
node, the emission of flashes of light in synchrony with the sound. 

The flash is emitted at the end of the sudden collapse the bubble suffers each acoustic cycle. This collapse takes 
around 4 /is, reduces the radius of the bubble from around 45 /im by a factor 10"^, and makes its surface to reach 
supersonic velocities. 

The violent deceleration of the bubble at the minimum radius is also accompanied by the emission of an outgoing 
acoustic pulse. After that, the bubble stays dead waiting for the next cycle. The posterior expansion occurs on 
hydrodynamic time scales, during the rarefaction half-cycle of the pressure (with some inertia making the bubble to 
reach its maximum size when the external sound field has already turned compressive). 

Each flash of light contains about one million of visible photons, and is approximately spherically symmetric. Its 
duration (less than 50 ps) is a hundred times shorter than the shortest (visible) lifetime of an excited state of the 
hydrogen atom. The energy of the photons ranges up to 6.5 eV (higher frequency photons cannot propagate through 
water), and the power of the flash can reach 100 mW. 

If the light were emitted from a region of atomic dimensions, a comparison of the flash energy with the average 
acoustic energy delivered to an atom of the fluid by the sound wave leads to conclude that a concentration of energy 
by twelve orders of magnitude should have occurred. 

This phenomenon is visible to the naked eye in a darkened room as a starlike light. 

Even though the hydrodynamical description of the collapse of the bubble and its posterior expansion is well 
understood |q], the mechanism through which part (about 0.01%) of the energy supplied by the sound is emitted as 
a flash of light is unknown and appears to be very complex ||5[] . 

Nevertheless, Schwinger suggested that the Casimir effect might be the underlying physics behind sonoluminescence, 
in the sense that the difference in the (static) electromagnetic zero point energy due to the change of the bubble's 
radius would be the available energy to be emitted as photons at the end of the bubble's collapse. In spite of the 
simplicity of this proposal, there has been no agreement about how to evaluate this change in the Casimir energy, 
and different approaches have led to controversial conclusions ||2l|-- |24| , [l7| , p5| --p8[ . 

In particular, the presence of singularities renders the physical interpretation of the energy finite parts difficult. 
However, this inconveniences may have their origin in the fact that the models usually employed in describing dielectric 
media mostly do not incorporate a realistic frequency dependent dispersion relation, then leading to an inadequate 
ultraviolet behavior. 

In a recent paper [ pi| , a nonmagnetic dielectric ball with a frequency dependent permittivity (a high energy 
approximation to the Drudc model), has been considered. It has been shown that a very simple pole structure results 
for the corresponding ('-function, and only a volume energy counterterm (to be absorbed in the mass density of the 
material) is needed to render the Casimir energy finite. Neither surface nor curvature counterterms are necessary. 

With the ultraviolet behavior under control, it makes sense to analyze the finite parts of the Casimir energy for 
realistic media. In this context, the analysis of the simple model under consideration is a step in the direction of 
incorporating finite frequency contributions. 

14 



For a spherical bubble of gas surrounded by water, we can take ni = 1 and n2 = 4/3. In this case, the pressure 
(Eq. (|66|)) reduces to 

Pia) /ny' _ -5 , 49 /an\-\ ^^/any'^ o.. .n-3 



, O ~ -2 X lO"-', (71) 

nn \cj 2167r2 54007r \ c ' ^ " ' ' ^ ^ 

approximately a negative constant if {afl/c) 3> 1, while the difference of Casimir energies (Eq. (|69|)) is 

Ecas.{ao)-Ecas.{a) rn\ ''^^^Vo-V)^2x IQ-^ (y^ - V) . (72) 

nil \c J 2lD7r^ 

Let us now consider an initial radius qq ~ 45 fim, and a final one a = ao/10. Then Vq ~ ^ = K)(l ~ 10"'^) ~ Vq. 

Firstly, we will estimate the difference of Casimir energies by equating it to the emitted energy. Assuming that the 
flash has one miUion photons with an average energy of 5 eV, we get (Ecas. (oo) — Ecas. (a)) ~ 5 10^ eV. Equation ( |7^ ) 
then gives (ao fl/c) ~ 608, which justifies the asymptotic expansion we have employed. The frequency cut-off turns 
out to be fi ~ 4 10^^ 1/s, equivalent to a (visible) energy of around 2.6 eV. Notice that the refraction index of water 
becomes essentially 1 at frequencies of the order 10^^ 1/s (see Ref. Q, page 291). The cut-off found corresponds to 
an electromagnetic pressure P ~ —2 10~^atm, of a much smaller magnitude than the acoustic pressure on the bubble 
(~ latm) 1). 

On the other hand, if we take instead as cut-off the frequency above which there is no propagation of photons in the 
water, fi ~ 10^® 1/s (corresponding to an energy of 6.5 eV, with oq fi/c ~ 1490), we get for the difference of Casimir 
energies (£'cas.(ao) — Ecas.{a)) — 1.8 10^ eV . The corresponding electromagnetic pressure is P ~ —7.5 lO^^'aim. 

Although obtained in the framework of a simplified model which ignores the complicated refraction index's de- 
pendence on the frequency, these results support Schwinger's proposal about the role the Casimir energy plays in 
sonoluminescence: It can behave as a reservoir of energy for the flash emission, which is feeded during the expansion 
of the bubble. 

VI. SUMMARY AND CONCLUSIONS 

In this paper we have considered a simple model of dielectric media, for which the permittivity and permeability 
are taken as constants up to a common cut-off frequency 17, above which they take the values corresponding to the 
vacuum. This assumption reflects itself in frequency dependent boundary conditions for the electromagnetic field at 
the interphase between dielectric materials, which becomes transparent for frequencies greater than the cut-off. 

For simplicity, we have limited our attention to nonmagnetic media, and studied the Casimir energy of a spherical 
dielectric of radius a and refraction index rii immersed in a second material of index n2, the whole contained in a 
large conducting concentric sphere of radius R. 

In this context, the (divergent) contribution of the frequencies higher than D, can be subtracted out by simply 
shifting the reference energy level. Indeed, it is independent of the (low frequency) refraction indices of both media, 
and also of the radius a of the internal sphere. 

On the other hand, the contribution of the eigenfrequencies lower than fJ reduces just to two finite sums, for 
transversal electric (TE) and transversal magnetic (TM) modes respectively. For each angular momentum I = z/— 1/2, 
these sums have been represented as differences of Cj^(s,2:) and (^J^^ {s,x), the incomplete ^-functions of the model 
(introduced in ||33| for the case of a scalar field with a frequency dependent boundary condition). 

In Appendix W we have derived the expression of the function A^^(z), whose roots determines the TE eigenfre- 
quencies for this configuration of material media. We have also proved that those among its zeroes lying in the open 
right half 2-plane are all real and simple. The same is true for A]^^^ (z). 

This fact allowed us to represent the incomplete ^-functions as integrals on the complex plane, employing the Cauchy 
theorem. Finally, the uniform asymptotic Debye expansion for Bcssel functions allowed for a systematic development 
of these C-functions, which facilitates the necessary analytic extensions. 

In Section O we have retained as many terms in this approximation as necessary to isolate the singular pieces of the 
incomplete C- functions at s = — 1. This has proved sufficient to evaluate the bulk and the first finite size contributions 
to the Casimir energy. 



Since the lowest positive zero of A^ (z) {A^ (z)) is a growing function of v, in Section HI we have determined 



Uq^ {vq'^), the maximum value of v for which there are eigenfrequencies smaller than or equal to fl. Starting from 
the analytic extension of the incomplete ^-functions to s ~ 0, we have been able to show that i^q^ and i^q^ are linear 
functions of (r7,2 R fl/c) (with corrections depending on lower non-integer powers of this parameter, induced by the 
presence of a cut-off), which are independent of ni and a. 
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In Section IV we have shown that, for {aVl/c) ^ 1, the dominant contributions are volume terms (see Eq. (p5|)). 
There is a piece proportional to the bubble's volume, whose sign is determined by the difference {n\ — rv^), plus a 
positive term proportional to the volume of the accesible space, which depends neither on ni nor on a. 

The second order in the Debye expansion produces the first finite size corrections to the Casimir energy, which are 
surface contributions (see Eqs. ( |60| ) and (|6l|)). There is a negative term proportional to the bubble's surface, coming 
exclusively from the TM modes. There are also surface terms corresponding to the external boundary, due to the TE 
and TM modes, which differ in sign and cancel out when added. 

Finally, there are also corrections proportional to non-integer powers of (n2 R Vt/c) (the power | and lower) induced 
by the relation between v^^ [vq^^] and R. These corrections, however, have no consequences in the R -^ oo limit. 

These results are analyzed in Section M. 

Firstly, we have considered the difference between the Casimir energy so evaluated and the one corresponding to 
the second medium filling completely the interior of the external sphere. As given in Eq. (^2|), this difference is 
independent of R. Its behavior with respect to a depends on the values of ni and 71,2 but, as remarked above, its 
values for different radius refer to different amounts of material media. 

In order to determine the force acting on the interphase between dielectrics, one should rather impose the conser- 
vation of the number of particles in each medium. This condition leads to a variation of the refraction indices with 
the bubble's radius a. 

In this way, we arrived at an expression for the electromagnetic vacuum pressure on the bubble, P{a), as a function 
of a (see Eq. (p9)). In this expression we could safely take the R ^ 00 limit to obtain the first terms of an expansion 
in powers of (afl/c)~^ ^ 1. P{a) so constructed vanishes for n2 — ni, and has a negative derivative with respect to 
n2. Therefore, it is negative for 712 > rii, and tends to compress the bubble. 

This pressure can be integrated to get the variation of the Casimir energy with respect to the bubble's volume, for 
given amounts of material media (see Eq. (|69|)). 

When considering models of dielectrics with constant refraction indices, the presence of divergencies makes it 
difficult to give a physical interpretation to the finite part of the vacuum energy, which cannot be isolated from the 
singular one. But when a realistic ultraviolet behavior is assumed, the singularities can be removed with a single 



volume energy counterterm |31|, to be absorbed in the mass density of the material. Neither surface nor curvature 
conterterms are needed to render the Casimir energy finite. With the ultraviolet behavior under control, one can 
worry about the finite frequency contributions in realistic models. 

The present paper, where low frequency refraction indices are modeled as constants up to the cut-off J7, can be 
considered as a step in this direction. 

Finally, we have applied the expressions found to a situation of interest for the phenomenon of sonoluminescence. 
Our results support Schwinger's proposal about the role the Casimir energy plays in the transduction of sound into 
light. 

Indeed, for the case of a spherical bubble of gas surrounded by water we can assume rii = 1 and n2 = 4/3. For a 
typical sonoluminescing bubble, the ambient radius is a ~ 4.5 fim (one tenth of its maximum radius). 

If we, moreover, estimate the difference in vacuum energies as the energy of a flash of light, we get from Eq. 
rif) an approximately negative constant electromagnetic pressure (favoring the collapse of the bubble, although of a 
magnitude much less than the acoustic pressure). Under these conditions, the cut-off f2 turns out to be in the region 
of the visible spectrum, large enough to justify the approach followed in this paper, and not very far from the region 
where the refraction index of water becomes essentially 1 . 

On the other side, if the frequency cut-off is imposed by hand where the propagation of light in water is no longer 
possible, then the change in the vacuum energy due to the collapse is about forty times the energy typically emitted 
in each flash. 

These results, obtained in the framework of a realistic dielectric model (which otherwise ignores the complicated 
refraction index's dependence on the frequency), clearly seems to support Schwinger's ideas about the role the Casimir 
energy plays in sonoluminescence: It grows with the bubble's volume by an amount comparable with the flash energy, 
which is therefore available to be emitted as light at the end of the collapse. 

Of course, this does not explain why the flash is emitted in such a short time at the end of the sudden collapse 
of the bubble. One could speculate about the formation of an excited electromagnetic field state, which would be 
induced to decay through some mechanism related to the strong deceleration stopping the bubble at its minimum 
radius. 
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APPENDIX A: EIGENFREQUENCIES OF THE TE MODES 

In this Appendix we will study the solutions of Eq. (P), subject to the boundary conditions for the TE modes, Eqs. 



(|6|). We will derive the expression of the function AfJ^^,^{z) in Eq. (|9|) (whose roots determine the eigenfrequencies) , 
and show that all its zeroes lying in the open right half z-plane are real and simple, a condition allowing for the 
integral representation in Eq. (|l5|). In particular, this implies that the only degeneracy of the eigenfrequencies is 
(21 + 1), due to the spherical symmetry of the problem. 

To this end, it is convenient to define 

s = sir)^l /!"' "^° (Al) 

[ f^2r + ai - a2, r > a 

with ai,2 — /ii,2 o- Then, expressing (p{s) = rfi{r) in terms of the new variable, taking into account that -^ = 
f ^ = Mfe II' ^ith fc = 1 (fc = 2) for r < a (r > a), and calling 

e(s) = ei9(ai - s) + 628(5- ai), ^(s) == /ii8(ai - s) + ^20(s - ai), (A2) 

we get, from Eq. (ph, the differential equation (with discontinuous coefficients) 






A ^{s) ^ ^ <! -^ - , . ^ ;:: — ^ ^is) = -^^(5), (A3) 



for s y^ ai- Here, / = 1, 2, . . . 

Moreover, Eq. (p3h results in the continuity conditions 



1^9(5 = a{) = ip{s = a^ ), and (p'{s = a^) = Lp'{s = ai ). (A4) 



So, we are looking for solutions of Eq. (A3) with a continuos first derivative, Lp{s) e C^(R+). 



In Subsection Al we will show that Eq. (A3) for this kind of functions, complemented with adequate boundary 
conditions, defines a self-adjoint operator. This excludes the pos sibility of non-real eigenvalues —lo^Jc^. The function 



A^ (z) is obtained in Subsection |A 2| . Finally, in Subsection A3 we will show that the nonvanishing zeroes of A^^ (z) 
are simple. 

1. Self-adjointness 



Let us consider the operator Ci defined as the differential operator Li in the left hand side of Eq. (A3), with a 
domain restricted to V{Ci) = C^[0, sq] C L^(R+, (e/^) (s) ds), where sq — ^2-R+ai — 02 (with R > a). Here C^[0, sq] 
is the space of functions with continuous derivatives of all orders and identically vanishing on some neighborhood of 
and Sq. Clearly, Ci is symmetric on 'D{Ci), 

(¥'l,'CiV?2)L2(R+,(e/^)(s) ds)= ('^iV'2,<^l)L2(R+,(e/^)(s) ds)" (A5) 

It is straightforward to show that its adjoint |^, C], is defined on the subspace 'D^Cj) C L^(R+, (e//x) (s) ds) 
containing those functions ip{s) with an absolutely continuous first derivative, and such that ip"{s) — Vl{s)^l){s) € 
L^([0, (5], ds), for 6 > (without requiring any further boundary condition). Moreover, for ijj G I?(£J), the action of 
CJ reduces to the application of the differential operator Ci . 

In order to determine the deficiency indices of Ci (defined as n±{Ci) = diTaKer{Cl + 1), see pTJ) one must look for 
the linearly independent solutions of CJ ?/'(s) = ±tip{s) in VlC]). 

Notice that the second derivatives of such functions are continuous for s 7^ ai. Moreover, if V'(s) is a solution of 
CJ ^(s) = +iiIj{s), then its complex conjugate ip{s)* is a solution oi CJ ^(s)* = —iiIj{s)*. This implies that Ci has 
equal deficiency indices, n_{Ci) = n_(_(£;), thus admitting self-adjoint extensions |37| ]. 

In fact, it can be seen that the equation 

17 



£J ip{s) = i^{s) 



(A6) 



has a unique (up to a constant factor) solution in !?(£[). Moreover, it vanishes at the origin. Therefore (see |p7|), 
there exists a one parameter family of self-adjoint extensions of Ci, which are in a one to one correspondence with the 
unitary maps from Ker{CJ — i) onto Ker{CJ + %), given by ?A(s) — > aip{s)*, where ip{s) is the solution of Eq. (|Aq ), 
and a G C with |a| — 1. 

Each essentially self-adjoint extension C^' is defined on a domain given by p^] 

P(/:[")) = {ct>{s) = ^(s) -I- /? [^(s) -H aV(s)1 ; ¥^(5) e Co°°[0, so], /3 G C}, (A7) 

with /:["^ acting on 0(s) e I?(/:|"^) C V{c\'''^) as 

4"V(s) = 4<^(s) = A<^(s) + z/3 (^(s) - a^(s)*) . (A8) 

In particular, notice that 0(0) = 0. 

Each essentially self-adjoint extension of Ci can also be characterized by the homogeneous boundary condition the 



functions in 'D{Ci') satisfy at s = sq- In fact, for all /? 7^ 0, Eq. (|A7|) imphes that 



'/>'(so) +c{a)(j){so) =0, 
with c{a) e RlJ{oo} (condition also satisfied for (3 — 0). 



(A9) 



2. The eigenfrequencies 

As seen before, we should impose on the fields a local homogeneous boundary condition at s = sq (i.e. r — R > a) 
This determines the functions to be included in the domain of the relevant operator, Eq. (A7). 

We choose to enclose the system within a large conducting sphere of radius R, obtaining the Dirichlct condition at 
s = Sq for the functions in the domain of an essentially self-adjoint extension of Ci which we call C 



(D). 



Eg,, 



cp{s) 



r=R 



= 0, V/>1. 



(AlO) 



So, the eigenfunctions (jy^^ of C\ satisfy the differential equation (A3), Ci(j)uj{s) = —{lu/c)4}^{s) for s 7^ ai, and 
the boundary and continuity conditions 



</>c.(s)L^o = o, 0c.(s)L=,„ = 0, 

This reduces the problem to looking for functions with a continuos second derivative for s 7^ ai, which satisfy 



(All) 



dzf 



1(1 + 1) 



4>ujis) = 0, 



(A12) 



with zi — s(a;/c)-y/ei/^i, for s < ai, which are solutions of the same equation with Zi ^ Z2 = (s— ai-f 02) (w/c) \/ ^ij [i^ , 
for s > fli, and satisfy the boundary and continuity conditions stated in Eq. (All). Therefore, 



(l)^{s) = AiJi+i/2{zi) + Biyi+i/2{zi), for s < ai, 

(j)^{s) = A2Ji+i/2{z2) + B2yi + 1/2{Z2), for Oi < S < Sq, 



(A13) 



where J^i+i/2{z) = z ji{z) and yi+i/2{z) = z yi{z) are the Riccati - Bessel functions, ji{z) and yi{z) being the spherical 
Bessel functions. 
Since 



J+i 



r(20 



Ji+i/2{z) = YmTW) (^ + ^("')) ' ^'+1/2^") = -^ (^ + ^("')) 



(A14) 
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the condition (j){0) = implies that Bi = 0,VL 
To ensure that 4'{so) = we can take 



A2 = ^^(+1/2(^0), B2 = -Ji+i/2{zq), (A15) 



with zo = (so - ai + 02) (w/c) ^62/^^2 • 

Finally, the continuity conditions at s = ai give 



AlJl+l/2{zi) ^yi+l/2{zo)Jl+l/2{z2) - Jl+l/2{zo)y 1+1/2(^2), 



^1^1 + 1/2(^1) "v/t"^ [yt + l/2{zo)J!+i/2{z2) - J/+l/2(^0)X+l/2(^2)} , 

where zi,2 = a(cj/c)^ei, 2^1,2- 

So, defining z = acu/c, the eigenfrequencies are determined by the zeroes of the function 



(A16) 



ATE 
^1+ 



V2(^) = = Jl+l/2{zi) [yi+l/2{zo)Jl + l/2(z2)- Jl + l/2{zo)y'i + i/2{z2)j 
^Jl+l^Zl) {yi+l/2izo)Jl+l/2{z2)- Jl+l/2{zo)yi+l/2{z2)} , 



(A17) 



where £ , , 






<TE /_"! j„f„,.,^:„„„ „,i „;„^,T,r,^,if^,. ^f r'^' 



Notice that every zero of A^5l/2(-^) determines an eigenvector of £j . Indeed, the function (l)uj{s) constructed as 
above is in I'(£j ) and satisfies Eq. (A3), —{uj/cY being the corresponding eigenvalue. 



Consequently, all the zeroes of A^ ,2(2) are either real or purely imaginary, since the operator £j is essentially 



self-adjoint for every ^ = 1,2, 



3. The multiplicities 



The condition determining the eigenvalues can also be understood from the following point of view. First notice 
that, for a given A: G C, the differential equation 

ti(t){s]k) = -k^(l){s-k) (A18) 

has, for s € (0,ai) and for s g (ai,so), two linearly independent solutions in C°°. So, a given function which is a 
solution at one side of the point s — ai can be continued as a solution to the other side, to obtain a C^(0, sq) function 
with a piecewise continuous second derivative. 



Now, let us call tp{s; k) a sq uare integrable non-trivial solution of Eq. ( [A18| ) satisfying (^(0; fc) = (unique up to a 



constant factor, see Eq. ( |A14| )): 

^{s)=Ji+i/2{zi), forO<s<ai, (A19) 



with zi = s k^Jei/ jii, and 

if{s) = A(fc)Ji+i/2(z2) + B{k)yi+y2{z2), for fli < s < So, (A20) 



with Z2 — {s ~ ai -\- 02) k\Je2/ 1^2- The coefficients A{k) and B{k) are determined by the continuity conditions at 
s = ai, 

A{k)Ji+y2{z2) + B{k)yi+y2{z2) - Ji + l/2(^l), 

, (A21) 

^(fc)^/+l/2(^"2)+5(fc)X+l/2(^"2) = ySf ^/+l/2(-^l). 

where zi,2 = a fc^ei, 2^^1,2- 

Similarly, let xi^', k) and p{s; k) be C^(0, sq) functions satisfying Eq. ( |A18| ) for s ^ ai, and the conditions 
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x(so;fc) = o, x'(so;fc)y^o, 
p{so;k) ^ 0, p'{so;k) = 0. 
For s > ai they can be taken as 

X{s;k) = yi+l/2{zo)Ji+i/2{z2) - Jl + l/2{zo)yi + l/2{z2), 
pis; k) = y'i^^/2i^o)Jl + l/2{z2) - Jl+i/2{zo)yi+l/2{z2), 



(A22) 



(A23) 



where zo = (s - ai + a2) k^/e^JJI^. 

Since all these functions are C^(0, so) solutions of Eq. ( Al8| ), the Wronskian of any two of them is a (/c-dependent) 
constant (for s < ai and for s> ai, and therefore for < s < Sq), which vanishes if and only if the selected functions 
are linearly dependent. In particular, 



W[p{s; k), x(s; k)] ={p(s; k)x'{s; k) - p'{s; k)x{s; k)}\^ 



'£2 



-kJ^{W[Ji+i,2{zo),yi+i/2{zo)]) y^O, forfc^O, 

V M2 



(A24) 



since Ji+ii2{z2) and 3^/4-1/2(^2) ar6 linearly independent solutions of Eq. ( A18| ) for s > ai. 
Let us call 

'q{k) = W[ip, x] = fis; k)x'{s; k) - if'{s; k)x{s; k), 
a{k) = W[ip, p] = ip{s; k)p'{s; k) - ip'{s; k)p{s; k). 



(A25) 



Since there are only two linearly independent solutions of Eq. ( |Al§| ), (f{s] k) can be expressed as a linear combination 
of x(s; fc) and p{s; k). In fact, 



V{k)p{s; k) - (j{k)x{s; k) ={(p{s; k)x'{s; k) ~ ip' {s; k)x{s; k)} p{s; fc)- 
- W{s] k)p'{s; k) - ip'{s; k)p{s; k)} x(s; fc)= W[p, x] ^{s; k). 



(A26) 



Consequently, (p{s] k) and x(s; k) are proportional for a given k (and, therefore, are C^(0, sq) solutions of Eq. (AIS) 
satisfying the Dirichlet boundary condition at s = 0, sq) if and only if ry(fc) = 0. 
Then, if ?7(fco) = we have 



^{s;ko) e'D{Cf'), and Cf\{s;kQ)^-klLp{s;ka). 



(A27) 



,{D) 



Since CY^^ is essentially self-adjoint, then /cq G R, which implies that the zeroes of rj{k) are either real or purely 
imaginary. 

Moreover, since rjik) is independent of s, the Wronskian can be evaluated at s = ai obtaining 



Vik) = W[lp, x]{s = ai) == Lp{a^ ; /c)x'(a+; fc) - ip'{a^ ; /c)x(a]^; A:), 



(A28) 



which is proportional to lS,Jr^,^{ak), as can be easily verified (see Eqs. (A19), (A23) and (A17)) 



We will finally show that the non- vanishing zeroes of ri{k) are simple. To this end, first notice that (as a fun cti on o f 
s) {dip/dk){s; k) £ C^iO, so), and {dtp/ dk){s; k) ~ s'+i for s -> 0+. This can be easily shown from Eqs. ( |aT9| ), ( [A20| ) 
and ( [A21 ). Moreover, from Eq. ( A18| ) we get 



i-'^W'M. 



"'"^11!*-')--^^*^-) 



(A29) 



for s ^ ai. 

Similarly, {dx/dk){s; k) g C-^(0, Sq) and, from Eq. ( A23 ), one can show that {dx/dk){s\ fc) ^ for s ^ sq. 

Let us n ow su ppose that fco 7^ is a multiple zero of rj{k). Then, ?7(fco) = and ri'{ka) = 0. It follows from Eq. 
( [A2^ ) and ( |A2^ ) that 
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— (s; fco) = Kix{s; ko) + K2-^{s; fco), 

for some constants Ki and K2- Therefore, we also have that {dip/dk){s; fco) -^ for s -^ sq. 
Consequently, {dtp/dk){s; fco) G T^i^i )• Moreover, since fc§ G R, from Eq. ( A29| ) we can write 



-2fco II (pis;ko) ||l2(r+. 



(R+,(e/Ai)(s)ds)- 



(p(s;fco), 



£; + ^0 



^(.;fco) 



£; + ^0 



f{s;ka),-^{s;ko) 



L2(R+,(£/p)(s)ds) 

= 0. 



(A30) 



(A31) 



L2(R+,(e/p)(s)ds) 

But this is a contradiction, since (^(s; fco) ^ 0. Therefore, all the non-vanishing zeroes of ?7(fc) are simple. 

An entirely similar analysis can be carried out for the TM case (now imposing the Neumann boundary condition 
at r = i?), to conclude that the roots of Af/^^,^{z) in Eq. (pi] ) lying in the open right half z-plane are all real and 
simple. 



APPENDIX B: DEBYE EXPANSIONS 



For the TE modes, making use of the uniform asymptotic expansion p^] of the Bessel functions appearing in the 
expression of A'^^{z), Eq . (U7h , with z ^ tvt (for i^ ^ 1 and t fixed), and discarding terms vanishing exponentially 
for R ^ CO, we get Eqs. ([J3D and (G4) where 



D^Tkt)=] (Vl + niH^-VT 



+ n2H^ + \ 1 + 



7122 i?2 i2 



(Bl) 



and 



4°i(t) -i ^ ' 



2t\l + ni^t^ l + n2^fi a^ + n2^RH^ VTTm^l^ a/TT^^^TF 



ta^n2^R^ Ua'-n2^RH^) J I 



712^ R^t^ 



t 



8(l + ni2t2)2 (l + „22t2) 
t 



8 (a2 + ^22 i?2 f2)3 
- [-8 7l2^+ni^i2 (^n2*t^-3~lO?22'^f) -4:711^ {2 + 7722^ f +n2*t^)] 



8 {1 + ni^ Py {1 + 712'^ f^)^ 



[ni" 712^ t" (n2^ t^ - 4) -712^ (8 + 3 772^ t^) - 2 77l2 (4 + 14 772^ t^ + 5 772^ t^ 



(B2) 



(B3) 



Similarly, for the TM modes, the use of the uniform asymptotic Debye expansion of the Bessel functions leads to 

dlog AP^(7i/t) 



dt 



= Dt''{t) + 0{v-% 



where D^^ (t) is given in Eq. (|2l) in terms of the algebraic functions 



^TmW =7 [Vl+m2i2_v/l + r722i2 + ^1 



n2 



^R^t^ 



(B4) 



(B5) 



J-'tm V') 



n^^t 



n2^t 



ni^n2^t 



2(l+77i2t2) 2(1+7722^2) j^^2 _^ ^^2 + 77l2 7722 ^2 2t 

ni^n2^t n2^R^t 



^/TTVi^^/TT^^7F (7712 + 7722 + 77i2 7722 1^) 2 {0? + 7732 i?2 1^) ' 



(B6) 
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and 



^TM V') — , ,5 

8a4(l + Ha!«ili)^ 



(B7) 



+ (n^^ t (-4n24 + n^' ^ + n,^ t^ (l + na^ t^^ (8 + na^ fi) 

-2ni^n2^ {6 + lln2^t^)+2niH^ {l + n2H^) (S + 13 ns^ t^ + 8 na^ i*) + 

ni4 na^ ^2 (-3 + 6 ^2^ i^ + Un^^ t^))) / (8 (l + rii^ ^2)^ (i + ^^^ ^2^1 („^2 ^ ^^2 ^ ^^2 „^2 ^2)2^ _ 

(ni^ i (8 712^ t^ (2 + 7122 t^) + ni^ n2^ (-12-3 n2^ t^ + 42 712'' f^ + 17 712^ t*^) + 
m^ (t2 + 14 7^2^t^ + 16 7^2^t*+7^2*i^'')+ 2 711^ (-2 + 7^22f2 (-11 + 3 7^22 f^ + 21 772'' t'' + 5 772^ i^)))) / 

(8 (1 + 77l2 t2) ^ (1 + 772^ t^)' (774^ + 772^ + Tll^ 772^ t^)^^ . 

In eqs. (B5-B7) we have also discarded those contributions vanishing exponentially for i? — > oo, which come from 
those terms containing K^i "^^"* ) in A^*^(z). 
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